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Notation 


Angle  of  incidence 
Circulation 

Circulation  density  at  a  given  angle 
Angle  of  strut 
Separation  angle 


deg 

m*/sec 

ro*/sec/rad 

deg  or  radians 
deg  or  radians 


Centre  of  total  circulation 
position 

Vorticity 
Maximum  vorticity 


angular 

deg  or  radians 

rad/sec 

rad/sec 


.  Coefficients  in  empirical  equations  various 

Body  diameter  metres 

Length  of  body  metres 

Induced  velocity  at  radius  from 

vortex  centre  metres/sec 

Radius  from  body  axis  to  position  of 

probe  metres 

Local  vortex  core  radius  (radius 

for  maximum  induced  velocity)  metres 


Non-dimensional  rate  of  turn  or  turn 
parameter;  L/S 

Local  radius  of  body  metres 

Centre  of  local  vortex  -  radial 

position  metres 

Centre  of  total  circulation  - 

radial  position  metres 

Radius  from  the  centre  of  the  vortex  metres 

Non-dimensional  radius  from  the 
centre  of  the  vortex:  R  /r* 


Arm  radius;  radius  of  turn 


metres 


Free  stream  velocity  in  the  X  direction  metres/sec 


Velocity  induced  in  Y  direction 

Velocity  induced  in  Z  direction 

Longitudinal  distance  from  nose  of 
body:  positive  aft 

Non-dimensional  distance  from  nose 
to  body;  x/L 

axis  which  the  vortex  lies  along 


Transverse  Planes 


metres/sec 

metres/sec 

metres 


ANALYSIS  OF  EXPERIMENTS  TO  MEASURE 
BODY  VORTICES  IN  CURVED  FLOW 


By  B  Ward 

A  R  J  M  Lloyd 

1 .  INTRODUCTION 

Experiments  on  the  Rotating  Arm  in  the  Manoeuvring  Tank  at  ARE(H) 
provided  data  on  vortices  shed  from  a  body  of  revolution  in  curved 
flow.  The  experiments  in  support  of  SUBSIM  (Reference  1)  took  place 
during  1986  (Reference  2)  and  1987  (Reference  3).  Reference  4  also 
gives  an  extensive  discussion  of  these  experiments.  This  memorandum 
describes  a  detailed  analysis  of  the  experiment  results. 

2.  VORITICITY;  SOME  BACKGROUND  READING 

2.1.  Characteristics  of  Vortices  (See  Reference  4) 

A  vortex  is  shown  in  Figure  1  with  the  free  stream  velocity  u  in  the  X 
direction.  Velocities  v  and  w  are  induced  in  Y,Z  directions 
respectively.  At  a  radius  R^  =  /y*+Z*  metres  from  the  centre  of  the 

vortex  a  velocity  q  =  Yv*+w*  is  induced  (see  Figure  2a).  From 
classical  potential  theory  q  =  r/27rRy  metres/sec  which  leads  to 

infinite  Induced  velocity  and  shear  at  the  origin  as  shown  in 
Figure  2b. 

In  reality  viscous  effects  stop  these  large  velocities  from  occurring. 
It  is  assumed  that  the  central  area  of  the  vortex  rotates  more  or  less 
as  a  solid  core  with  the  induced  velocity  close  to  the  centre 
increasing  linearly  with  the  radial  distance  from  the  centre.  This  is 
simulated  by  the  modified  equation 


q=  ^ 

1  -  exp 

-  1.26  R’2 

2nR 

V 

V  J 

metres/second 


where 


and  r*  is  the  core  radius  of  the  viscous  vortex  (defined  as  the 
radius  at  which  the  induced  velocity  q  is  a  maximum) . 

Here  the  induced  velocity  is  reduced  to  zero  at  the  centre  but  reaches 
its  inviscld  value  at  large  distances  from  the  centre. 


1 


2 


2.2.  Body  Vortices 

If  the  body  is  at  a  aoderate  angle  of  incidence  and  turn  rate  the 
boundary  layer  on  the  leeward  surface  separates  as  shown  in  Figure  3. 
The  location  of  the  separation  line  may  be  quantified  in  terws  of  the 
'separation  angle'  0^.  Vorticity  shed  from  the  boundary  layer  is 

convected  away  and  coallesces  to  form  a  diffuse  pair  of  vortices  with 
cores  almost  parallel  to  the  body  axis.  The  vortices'  strength 
increases  towards  the  tall  of  the  body  as  more  vorticity  is  added. 
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3. 


THE  EXPERIMENT 


A  5.0  metre  long  model  body  of  revolution  was  attached  to  the  Rotating 
Arm  as  shown  in  Figure  4a. 

Table  1  gives  the  model  dimensions. 

Table  1 

MODEL  DIMENSIONS  (LENGTH.  5  m;  MAXIMUM  RADIUS,  0.294  m) 


X’ 

R  (m) 

0 

0 

0.05 

0.206 

0.10 

0.266 

0.15 

0.291 

0.20-0.65 

0.294 

0.70 

0.289 

0.75 

0.274 

0.80 

0.247 

0.85 

0.206 

0.90 

0.151 

0.95 

0.084 

1.00 

0 

A  Freestone  vorticity  probe  (Reference  5)  was  designed  and  built  with 
a  radius  of  10  mm.  The  probe  was  attached  to  a  strut  which  was  built 
to  provide  measurements  out  to  a  radius  of  1.5  x  body  diameter.  An 
inclinometer  was  used  to  set  up  the  strut  such  that  in  the  position 
referred  to  as  0  =  90  degrees  it  protruded  vertically  from  the  top  of 

the  body  and  at  0  =  0  degrees  it  was  horizontal  on  the  inboard  side  of 

the  model,  (see  Figure  4b). 

Three  different  strut-body  stations  were  used  in  the  test  conditions 
(a)  x'  =  0.925,  (b)  x'  =  0.85,  (c)  x'  =  0.7  where  x'  =  x/L  (see 
Figure  4b ) . 

Incidence  angle  a  was  set  with  the  nose  towards  the  centre  of  the  arm 
(see  Figure  4c). 

The  probe  was  traversed  from  close  to  the  body  surface  to  the  limits 
of  the  strut  where  the  vorticity  was  expected  to  be  zero.  This  was 

done  with  the  strut  at  angles  of  0,  10,  20.... 90  degrees.  The  arm 

was  rotated  counter-clockwise.  The  non-dimensional  rate  of  turn  or 
turn  parameter  r'  was  varied  by  changing  the  arm  radius  S.  About  130 
data  readings  were  obtained  for  each  test  condition. 
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4. 


TYPICAL  DATA  SET 


THIS  REGION 


Figure  5 

TYPICAL  VORTICITY  TRAVERSE 

Figure  5  shows  a  sketch  of  a  typical  vorticity  traverse  (plot  of 
vorticity  versus  radius  for  a  given  angle  0).  A  number  of  salient 
features  are  apparent; 

a.  The  peak  value  has  probably  been  missed  due  to  the  limited 
number  of  data  values. 

b.  The  physical  dimensions  of  the  probe  precluded  any 
measurements  close  to  the  local  body  svvface. 

c.  The  traverses  appear  to  come  down  to  a  'plateau'  as  the 
radius  increases,  whereas  it  would  have  been  expected  that  they 
would  decrease  smoothly  to  zero  at  the  'edge'  of  the  vortical 
flow.  This  is  probably  associated  with  zero  errors  in  the 
Freestone  probe.  It  was  found  in  Reference  2  that  the  probe 

gave  small  levels  of  'vorticity'  in  uniform  flow  when  no  vorticity 
was  expected  and  that  these  zero  errors  were  functions  of  the 
flow  direction  at  the  probe.  A  simple  correction  was  devised 
based  on  the  results  obtained  at  the  outer  limits  of  the  strut, 
and  this  was  applied  to  each  traverse,  forcing  the  vorticity  at 
r  =  1.5  X  D  to  be  zero.  However,  the  flow  direction  changes  as 
the  probe  is  moved  towards  the  body  and  this  will  affect  the  zero 
error  correction.  No  simple  method  of  quantifying  this  effect 
has  been  devised. 

To  overcome  these  shortcomings  the  data  were  faired  to  a  more  expected 
form  and  ensuring  smooth  variations  as  each  of  the  parameters  was 
varied.  This  Involved; 
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a.  Inserting  estiaated  additional  points  close  to  the  local 
radius  of  the  body. 

b.  Interpolating  extra  points  throughout  the  curve  including, 
at  times,  a  point  for  peak  vortlcity. 

c.  Bringing  traverses  down  to  zero  at  the  apparent  edp'^  of  the 
vortical  flow. 

d.  Ensuring  that  the  estimated  peak  vorticlty  increased 
smoothly  with  strut  angle. 

5.  DATA  ANALYSIS  PROGRAM 

5.1.  Introduction 

A  FORTRAN  program  entitled  INTER?  was  developed  on  the  ARE(H)  PRIME 
9955  Computer  to  fair  the  data.  The  fairing  involved  human  judgement 
in  altering  and  inserting  data  points  depending  on  the  trend  of  the 
curve  and  the  expected  behaviour  of  the  traverses. 

The  program  used  the  GINO  graphics  package  which  enabled  plots  of 
vorticlty  traverses  to  be  displayed  on  the  VDU,  analysed,  subsequently 
altered  and  analysed  further. 

5.2.  Interpolating  Extra  Data  Points 

Each  vortlcity  traverse  (9  =  10  degrees  to  90  degrees)  was  analysed. 
When  it  was  considered  appropriate  cubic  polynomials  were  calculated 
to  fit  all  sets  of  four  adjacent  points  (1.2, 3, 4),  (2, 3, 4, 5)  etc,  for 
the  particular  angle  under  examination.  This  involved  taking  the  four 
points  (Figure  6)  deriving  a  system  of  equations  and  then  solving 
them.  The  resulting  cubic  polynomial  C(r)  =  a^r^  +  Ugr*  +  a^r  +  a^ 

was  used  to  calculate  an  additional  point  between  the  middle  two 
points ,  (Figure  7 ) . 


Figures  6  and  7 
INTERPOLATION  OF  EXTRA  POINT 
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Doing  so  throughout  the  curve,  N-3  extra  points  were  interpolated 
(where  N  is  the  total  nunber  of  original  data  points).  These  extra 
points  ensured  a  more  accurate  numerical  integration  and  in  certain 
cases  a  more  reasonable  peak  value  for  vorticity  was  achieved. 

5.3.  Additional  Points  Close  to  the  Body  Surface 

Figures  8  and  9  show  how  a  traverse  for  a  strut  angle  of,  say 
30  degrees,  was  extrapolated  to  the  body  surface.  In  this  example  the 
trend  of  the  curve  shown  implies  that  the  vorticity  is  negative  at  the 
body  surface. 


Figures  8  and  9 

EXTRAPOLATING  EXTRA  POINTS  CLOSE  TO  BODY  SURFACE 
FOR  STRUT  ANGLE  9a  30  DEGREES 


Figures  10  and  11  show  how  a  traverse  for  a  high  strut  angle  of,  say 
90  degrees,  was  extrapolated  to  the  body  surface.  At  such  a  strut 
position  it  was  assumed  that  the  vorticity  would  peak  at  the  body 
surface.  Hence  the  trend  of  the  curve  will  be  as  shown  in  Figure  11. 
The  extrapolated  maximum  vorticity  at  the  body  surface  at  angles  of 
80  degrees  and  90  degrees  was  subsequently  reassessed  by  the  method 
described  in  5.5. 


EXTRAPOUTING  EXTRA  POINTS  CLOSE  TO  BODY  SURFACE 
FOR  STRUT  ANGLE  9«  90  DEGREES 
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Elialnatlng  the  Zero  Plateau 


r 


Figure  15 

ADJUSTED  VALUES  OF  PEAK  VORTICITY 
(x’  =  0.925,  r’  =  0.3,  a  =  7,5  DEGREES) 

5- 5-  Increasing  the  Peaks  of  Vortlcity  Traverses 

A  regression  polynonial  of  the  for«  =  bQ+bj0+b20*+b30=>+b40*+b50s 

was  fitted  to  the  estimated  peak  vortlcity  data  for  each  test 
condition.  An  example  is  shown  in  Figure  14.  It  was  argued  that  the 
peak  vortlcity  should  Increase  monotonically  up  to  a  maximum  at  about 
0  =  90  degrees  where  the  boundary  layer  separates  from  the  body 
surface.  If  this  is  the  case  the  estimated  peak  values  at 
0  “  80  degrees  and  90  degrees  are  clearly  too  low.  Recalling  that 
these  values  were  estimates  based  on  extrapolation  (as  shown  in 
Figures  10  and  11)  it  was  considered  Justifiable  to  Increase  their 
values  as  shown  in  Figure  15. 

The  regression  curve  was  used  to  assess  whether  any  further 
adjustments  to  the  estimated  peak  levels  were  required  at  other  values 
of  0.  In  Figure  15  an  example  of  such  an  adjustment  can  be  seen  for 
0  >  40  degrees. 

In  certain  cases  these  changes  demanded  reassessments  of  earlier 
decisions.  In  particular  increasing  the  peak  vorticity  levels  at  high 
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Figure  16 

FAIRED  VORTICITY  TRAVERSES  -  AN  ARTIST  S  IMPRESSION 


BEFORE 


AFTER 


RADIUS  MM  X  10* 

Figure  18c 

VORTICITY  TRAVERSES  -  BEFORE  AND  AFTER  FAIRING 
’  -  0.925.  r’  ° 


values  of  9  often  aade  reassessaent  of  the  position  of  the  edge  of  the 
vortical  flow  necessary. 

An  artist's  lapression  of  the  fairing  of  vortlclty  traverses  can  be 
seen  in  Figures  16  and  17.  For  further  exawples  of  the  effects  of 
fairing  on  traverses  see  Figures  18a  to  18c. 

5.6.  Circulation 

5.6.1.  Local  Circulation  Density 

In  Figure  19  the  circulation  residing  within  a  sector  defined  by  the 
saall  angle  dO  is 
00 

r^dO  =  /  C  r  dr  d0  m^/sec 
o 

and  the  local  circulation  density  is 
00 

Tg  =  /  C  r  dr  a^/sec/rad 
o 


This  local  circulation  density  was  calculated,  like  all  other 
integrals  in  the  analysis,  using  Trapezoidal  rule. 

Both  negative  and  positive  circulation  densities  were  calculated. 

As  a  zero  f-value  was  not  always  present  in  data  where  traverses  cut 
the  r-axis  on  the  inboard  side  of  the  curve  the  method  of  ' regula 
falsi'  was  used  to  calculate  the  intersecting  point  as  shown  in 
Figure  20. 

"2  “  ‘‘1  ■  -  Co 


REGULA  FALSI 
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CIRCULATION  DENISITY  m^/  SEC /RAD 


A  graph  of  the  circulation  density  against  angle  was  plotted  for  each 
experlaent  condition  and  an  exaaple  is  shown  in  Figure  21. 


CIRCULATION  DENSITY  AGAINST  STRUT  ANGLE 
(x'  =  0.925,  r'  -  0.3.  a  =  7.5  DEGREES) 


5.6.2.  Total  Circulation 

The  total  circulation  was  estimated  from 
71/* 


r  =  /  Fg  d0 
o 

Results  for  all  test 


m*/sec 

conditions  at  x' 


=  0.925  are  shown  in  Figure  22. 
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TURN  PARAMtTER  r* 


Figure  22 

TOTAL  CIRCULATION  AGAINST  TURN  PARAMETER 
FOR  VARYING  VALUES  OF  INCIDENCE  ANGLE 

5.7.  Vortex  Centres 

5.7.1.  Local  Vortex  Centre  -  Centre  of  Circulation  Density 

The  radial  location  of  the  centre  of  the  circulation  density  at  a 
given  angle  9  is; 

00 

/  f  r*  dr 

R.  »  - -  «etres 

1  00 

/  f  r  dr 

o 


(See  Figures  19  and  23) 


Values  were  calculated  for  all  test  conditions  and  an  example  is  shown 
in  Figure  24. 

5.7.2.  Total  Vortex  Centre  -  Centre  of  Total  Circulation 

The  centre  of  the  total  circulation  was  calculated  using  the  following 
integrals ; 


00  TT/  2 

/  J  C  1*“  0  (i®  dr 

(R6)i>tr  =  - — -  metres  .  radians 

OV  71/2  00 

/  /  C  r  dr  do 

o  o 

rr/a  «> 

/  /  C  dr  d0 

R„,,  =  ^ ; — 2 -  metres  (centre  of  circulation-radial 

DV  77/2  00 

/  /  c  r  dr  d9 

o  o 

(R0)bV 

=  -  radians  (centre  of  circulation-angular  position) 

•^BV 

Results  for  all  test  conditions  at  x'  =  0.925  are  shown  in  Figures  25 
and  26.  (For  convenience  9gy  is  plotted  in  degrees). 
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THETA  DEG 

Figure  24 


LOCAL  VORTEX  CENTRE  AGAINST  STRUT  ANGLE 
(X'  =  0.925.  r'  =  0.2,  a  =  7.5  DEGREES) 


5.8. 


Vortex  Core  Radius 


According  to  Reference  6  the  peak  vorticity  and  the  vortex  core  radius 
are  related  by  the  formula 


1.26  r 


‘MAX 


n  rj 


rad/sec 


Therefore 


r*  = 


n.2e  r 


e 


n  c 


MAX 


metres 


Figure  27  shows  core  radii  calculated  in  this  way  for  a  particular 
test  condition. 


Figure  25 

CENTRE  OF  CIRCULATION  -  RADIAL  POSITION  AGAINST 
TURN  PARAMETER  FOR  VARYING  VALUES  OF  INCIDENCE  ANGLE 
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THETA  DEG 


Figure  27 

VORTEX  CORE  RADIUS  AGAINST  STRUT  ANGLE 
(x'  =  0.925,  r'  =  0.3,  a  =  2.5  DEGREES) 


6. 


EMPIRICAL  EQUATIONS  TO  REPRESENT  THE  RESULTS 


6.1.  Circulation  Density 

Plots  of  circulation  density  against  strut  angle  for  varying  test 
conditions  were  examined  (see  Figures  21  and  28).  The  trend  of  the 
plots  suggested  a  sine  curve  would  be  appropriate  with  varying 
amplitude  determined  by  the  parameters  r'.  a  and  x'.  The  following 
equation  was  determined. 


C2r ' *  +  Cjr '  +  Cq  + 


d2r'*  +  djT'  -  do|a 


sin 


fg  0=-  ^  -  fl0 


•0  "  gi^ 


(x'^0.65)  m^/sec/rad 


Tq  =  0 


(x'<0.65) 


(q ,  a  in  degrees) 


The  coefficients  of  the  above  equation  are  noted  in  the  following 
table . 


•=2 

Cl 

Cq 

*^2 

^1 

*^0 

^3 

^2 

go 

gl 

-15.0 

9.3 

-1.4 

2.0 

00 

o 

1 

0.208 

0.00012 

-0.037 

4.23 

-2.398 

3.66 

NOTE: 

Fg  =  0  for  all  values  of  0  at  x’  =  0.65  where  circulation  becomes  zero. 
The  quality  of  the  fit  can  be  seen  in  the  example  shown  in  Figure  28, 

6.2.  Total  Circulation 

For  x'  =  0.925.  plots  were  made  of  the  circulation  against  turn  parameter 
for  each  angle  of  incidence  as  shown  in  Figure  22. 

The  values  appeared  very  consistent  for  2.5  degrees,  5.0  degrees  and 
7.5  degrees  angles  of  incidence.  As  data  for  0  degrees  and  10  degrees 
had  come  from  the  1987  experiment  when  the  probe  was  known  to  have 
been  damaged,  it  was  decided  to  base  any  further  mathematical  models 
largely  on  the  trends  of  the  three  middle  angles,  but  still  taking 
0  degrees  and  10  degrees  into  account. 

A  regression  polynomial  of  degree  5  with  respect  to  r'  fitted  the  data 
for  a  >  2.5  degrees,  5.0  degrees,  7.5  degrees  very  well  and  the 
required  behaviour  as  r'-  0  was  satisfied.  The  constant  term  of  the 
polynomial  was  varied  according  to  a. 
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Figure  28 

CIRCULATION  DENSITY  AGAINST  STRUT  ANGLE 
(x'  =  0.925,  r’  =  0.4,  g  =  5  DEGREES) 


The  equation  for  x'  =  0.925  is  as  follows. 


^(x'=0.925)  ■  ^5  h^r'^^  hg  r'’+  hg  r'%  hj  r'+  jgO' 

HI* /sec  (a  in  degrees) 

A  more  general  equation  valid  for  all  x'  is  as  follows. 


r  = 


*"(x'=0.925)  (x’>0.65) 


■*/sec 


0 


(x'<0.65) 


The  coefficients  of  the  above  equations  are  noted  fh  the  following 
table. 


“4 

^2 

hi 

m 

D 

D 

^0 

258.78 

-279.81 

99.38 

-11.81 

0.48 

0.004 

0.018 

3.67 

-2.39 

The  quality  of  the  fit  of  the  equation  at  x'  =  0.925  can  be  seen  from 
Figure  22. 

6.3.  Local  Vortex  Centres 

Plotting  the  local  vortex  centres  radial  position  against  strut  angle 
for  varying  conditions  it  was  Immediately  apparent  (see  Figure  24) 
that  a  linear  equation  was  appropriate. 

The  following  equation  was  derived. 


R 


1 


+  ((Po+Pl®  (qQ^-q^alr '  )0) .  (Sq*s^x' ) 

+R  (x'>0.65) 

R  (x'<0.65) 


metres 


(9  and  o  in  degrees) 

The  coefficients  of  the  above  equation  are  noted  in  the  following 
table . 


'o 

ll 

■"o 

mi 

"o 

"l 

30.864 

15.97 

212.89 

-  12.4 

-  0.00233 

0.00359 

Po 

Pi 

^0 

^1 

So 

Si 

-  0.344 

-  0.165 

-  1.857 

0 . 0397 

-  0.00229 

0 . 00355 

At  x'  =  0.65  where  the  circulation  becomes  zero  the  value  of  will 
equal  the  local  body  radius  for  all  values  of  9. 

The  quality  of  fit  (see  Figure  24)  of  the  equation  was  consistently 
good  for  all  test  conditions  with  x’  =  0.925. 

6.4 .  Total  Vortex  Centre 

The  radial  locations  of  the  centre  of  circulation  (for  x'  =  0.925) 
were  plotted  against  turn  parameter  for  each  angle  of  incidence 
(Figure  25). 
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A  regression  polynoalal  of  degree  3  was  found  to  fit  data  for 
a  =  2.5  degrees.  5.0  degrees,  7.5  degrees  with  x'  =  0.925  very  well. 
The  following  equation  was  developed. 


(x- 


=0.925) 


tg  r'*  +  tg  r'*  +  tj  r'  +  VgO*  +  VjO  +  0.118 
(a  In  degrees) 


metres 


and  nore  generally 


(exp(ZgX'®  +  ZgX'*  +  ZjX'  +  Zq)  +  Wg) 


.  (tg  r '  *  +  tg  r'  *  +  tj  r ' 

+  Vg  O*  +  Vj  a)  g  jjg 

+  R 

(x'>0.65) 

R 

(x'<0.65) 

metres 


(a  in  degrees) 

The  coefficients  of  the  above  equations  are  noted  in  the  following 
table . 


^3 

^^2 

D 

'^2 

''l 

-  1.353 

0.4492 

0 . 00064 

0.00296 

"3 

^2 

^1 

^0 

Wo 

5.06 

-  8.0 

5.62 

-  1.67 

-  0.993 

Note  that  at  x'  >0.65  the  radial  position  is  equal  to  the  local 
radius  of  the  body. 

Likewise  from  plots  of  the  0  position  (for  x'  =  0.925)  against  turn 
parameter  for  each  angle  of  Incidence,  (see  Figure  26),  the  following 
linear  equation  was  developed. 

0gy  -  Ajr'  +  BgO*  +  BgO*  +  Bjtt  +  Bg  deg 

(X'  -  0.925) 

(a  in  degrees) 
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A  general  equation  was  developed. 


e 


BV 


cos 


V* 


+  FjX' 


+  1 


.  Ajr'  ^  Bga^  >  ^  B^a  .  B^j  (x'>o.65) 

(a  In  degrees) 

The  coefficients  of  the  above  equations  are  noted  in  the  following  table. 


83 

^2 

®0 

^2 

^0 

7.76 

0.024 

-0.162 

-2.25 

47.68 

6.724 

-9.546 

4.648 

6.5.  Vortex  Core  Radius 

By  calculating  the  core  radius  over  angles  10  degrees  to  90  degrees 
for  all  test  conditions  and  by  analysing  the  resulting  values  a  linear 
equation  appeared  appropriate. 

The  following  general  equation  was  developed. 

r,  =  X'  [Gq  +  HqO  +  J^ej  metres  (x’>0.65) 

(a, 9  in  degrees) 

The  coefficients  of  the  above  equation  are  noted  in  the  following 
table . 


% 

«0 

•^0 

0.0526 

0 . 00095 

-0 . 00046 

The  quality  of  the  fit  can  be  seen  from  Figure  27.  It  should  be  noted 
that  the  equation  did  not  fit  all  the  test  condition  data  as  well  as 
that  of  the  example  but  as  mentioned  previously  not  all  the  data  can 
be  considered  very  reliable.  The  above  equation  was  derived  from 
analysis  of  the  most  consistent  data.  Nevertheless  the  equation  does 
reflect  the  linear  trend  and  magnitude  of  most  of  the  results. 

7.  CONCLUSIONS 

Experiments  on  the  Rotating  Arm  in  the  Manoeuvring  Tank  at  ARE(H) 
provided  data  on  vortices  shed  from  a  body  of  revolution  in  curved 
flow.  The  data  contained  apparent  shortcomings  as  explained  in  this 
report.  Methods  as  described  above  were  used  to  fair  the  data  to  a 
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■ore  expected  for*.  The  resulting  faired  data  was  used  to  develop 
e«plrlcal  equations  for  local  circulation  density,  total  circulation 

density,  apparent  centre  of  total  circulation,  local  vortex  centre  and 
local  vortex  core  radius. 

These  expressions  will  be  Incorporated  In  the  SUBSIM  Computer  program 
In  due  course. 
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